Deep inelastic prompt photon production at HERA 
in the fc^-factorization approach 



S.P. Baranov", A.V. Lipatov^ N.P. Zotov^ 
January 27, 2010 



" P.N. Lebedev Physics Institute, 
119991 Moscow, Russia 

^ D. V. Skobeltsyn Institute of Nuclear Physics, 
M. V. Lomonosov Moscow State University, 
119991 Moscow, Russia 



Abstract 

We investigate the prompt photon production in deep inelastic scattering at HERA in the 
framework of fc-r-factorization QCD approach. Our study is based on the off-shell partonic 
amplitude eq* — )■ ejq, where the photon radiation from the leptons as well as from the 
quarks is taken into account. The unintegrated quark densities in a proton are determined 
using the Kimber-Martin-Ryskin prescription. The conservative error analisys is performed. 
We investigate both inclusive and jet associated prompt photon production rates. Our 
predictions are compared with the recent experimental data taken by the HI and ZEUS 
collaborations. We demonstrate that in the /cT-factorization approach the contribution from 
the quark radiation subprocess (QQ mechanism) is enhanced as compared to the leading- 
order coUinear approximation. 



PACS number(s): 12.38.-t, 13.85.-t 



1 Introduction 

The production of prompt photon^ in deep inelastic ep collisions at HERA is a subject of 
intense studies [1-3]. Such processes provide a direct probe of the hard subprocess dynamics, 
since the produced photons are largely insensitive to the effects of final-state hadronization. 

^Usually the photons are called "prompt" if they originate from hard interaction, not from the decay of 
final-state hadrons. 
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Studying the photon production in deep inelastic scattering (DIS) provides a test of per- 
turbative Quantum Chromodynamics (QCD) with two hard scales: the transverse energy of 
the emitted photon Ej^, and the exchanged photon virtuality Q^. Good understanding of 
the production dynamics is also important for searches of new particles decaying to photons 
at the LHC conditions. 

At the leading-order, the production of prompt photons in deep inelastic ep scattering is 
described by the parton-level process eg — )■ e'jq, where q represents a quark or an anti-quark, 
and e represents the beam lepton. The corresponding cross section a{ep — )■ e'jX) can be ob- 
tained by convoluting the parton-level cross section cr{eq — )■ e'jq) with the quark distribution 
functions in the proton. The observed final-state photon can be emitted by a quark or by a 
lepton (the so-called QQ and LL mechanisms, respectively). The interference contribution 
(LQ) is small. The QQ mechanism includes the direct photon radiation from the quark line 
and also the so-called fragmentation processes [4], where the produced quark forms a jet 
containing a photon with large fraction of the jet energy. This contribution involves poorly 
known quark-to-photon fragmentation functions fq-^^{z) [5]. However, the isolation criterion 
introduced in experimental analyseisubstaufally reduces the fragme„tat,o„ eompouent |1- 
3], 

Recently the HI and ZEUS collaborations have reported their data [1-3] on the deep 
inelastic production of prompt photons at HERA, both inclusive and in association with 
a hadronic jet. The presented data on the inclusive prompt photon production have been 
compared [2, 3] with the LO pQCD 0{a^) calculations [6]. A substantial (by a factor of about 
2) disagreement between the data and theory has been found [3] at low Q^. Also, the absolute 
size of the reported experimental cross sections turned out to exceed the predictions of Monte- 
Carlo generators herwig [7] and pythia [8] (by the factors of 7.9 and 2.3, respectively). 
Even after normalizing the total production rate, none of these programs was able to describe 
all kinematical dependences of the measured cross sections [1-3]. Comparison between the 
LO pQCD calculations [6] and the data in bins of r]'^ shows that the difference can be 
attributed [3] mainly to the underestimation of the QQ contribution. At present, the next- 
to-leading order (NLO) 0{a^as) predictions are only available [9] for the associated photon 
plus jet phase space selection, yet not for inclusive process. These predictions are higher than 
the results of LO calculations [6], especially at low Q^, but still underestimate the data [3]. 
It was concluded [1-3] that further theoretical investigations (including the evaluation of 
higher-order processes) are needed to understand the observed discrepancy. 

In the present paper to analyse the HI and ZEUS data [1-3] we use instead the stan- 
dard (collinear) approximation the so-called /cr-factorization [10] approach, which is based on 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) [11] or Ciafaloni-Catani-Fiorani-Marchesini (CCFM) 
[12] evolution equations. A detailed description of the fcT^-factorization approach can be 
found, for example, in the reviews [13]. Recently we have applied this approach to the inclu- 
sive and jet associated prompt photon photoproduction at HERA [14] and hadropro duct ion 
at Tevatron [15]. As it was demonstrated in the ZEUS paper [17] and also in the recent 
study [17] performed by the HI collaboration, the ^T-factorization predictions are in better 
agreement with the data than collinear NLO calculations. 

In the present study we extend our investigations [14] to the DIS region. We calculate 

^See also Section 2.4. 
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the relevant off-shell (i.e. /cy-dependent) amplitude Ai{eq* — )■ e'-fq), where the virtuality of 
the incoming quarks is properly taken into account. Our consideration covers both the LL 
and QQ production mechanisms. We neglect the LQ mechanim since it gives only about 3% 
contribution to the total cross section. Then, we use the obtained expressions to calculate 
the total and differential cross sections of the deep inelastic production of prompt photons 
(both inclusive and accompanied by a hadronic jet). The unintcgrated quark densities in 
a proton fq{x,li^, n^) are taken in the Kimber-Martin-Ryskin (KMR) form [18]. We make 
a systematic comparison of our predictions to the available HI and ZEUS data [1-3]. Our 
main goal is to study the capability of the /cT-factorization approach in describing the DIS 
data. Such calculations are performed for the first time. 

The outline of our paper is the following. In Section 2 we recall shortly the basic formulas 
of the /cT-factorization approach with a brief review of the calculation steps. In Section 3 
we present the numerical results of our calculations and a discussion. Section 4 contains 
our conclusions. The compact analytic expressions for the off-shell matrix elements of the 
eg* — > e'yq subprocess are given in the Appendix. 



2 Theoretical framework 



2.1 Off-shell amplitude of the eq* — > e'yq subprocess 



Let us denote the four-momenta of the incoming electron and proton and the outgoing 
electron, photon and final quark by p^, Pp, Pe', Pj and Pq, respectively. The initial quark and 
virtual photon have four-momenta k and q. 

Let us start with the LL production mechanism. There are two simple Feynman diagrams 
which describe this partonic subprocess at the leading order in ccem- The relevant matrix 
element reads 

MLLieq* ^ e^q) = Cqeh^^'^^L^^H^ (1) 

where tq is the fractional electric charge of the quark and e'^ is the produced photon polariza- 
tion four- vector. The leptonic and hadronic tensors are given by the following expressions: 



7^ , , „ NO „ o 7 +7 



U{pe), 



(2) 



H^^u{pq)Yu{k). (3) 

When we calculate the matrix element squared, the summation over the produced photon 
polarizations is carried with X^e'^e'' = —g^^i and the spin density matrix for all on-shell 
spinors is taken in the standard form u{p)u{p) = p + m. In the case of off-shell initial quarks 
the on-shell spin density matrix has to be replaced with a more complicated expression. To 
evaluate it, we "extend" the original diagram and consider the off-shell quark line as internal 
line in the extended diagram. The "extended" process looks hke follows: the initial on-shell 
quark with four-momentum p and mass ruq radiates a quantum (say, photon or gluon) and 
becomes an off-shell quark with four-momentum k. So, for the extended diagram squared 
we write: 



iMI'^Sp 



k -\- TUq / V _ . s k + rUq 

— -1 u[p)u[p)-i^ 



k^ — rriq 



j- 

k"^ — mi, ^ 



(4) 
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where T is the rest of the original matrix element which remains unchanged. The expression 
presented between T'^ and 7^ now plays the role of the off-shell quark spin density matrix. 
Using the on-shell condition u{p) u{p) = p+rUq and performing the Dirac algebra one obtains 
in the massless limit rUq — > 0: 

l-^l' - pp'^'^ (2^'^ - 4(P • k)k) %. (5) 

Now we use the Sudakov decomposition k = xp + kr and neglect the second term in the 
parentheses in (5) in the small- a; limit to arrive at 

i-^r-^r-^pT,. (6) 

(Essentially, we have neglected here the negative light-cone momentum fraction of the in- 
coming quark). The properly normalized off-shell spin density matrix is given by xp, while 
the factor 2/xk'^ has to be attributed to the quark distribution function (determining its 
leading behavior). With this normalization, we successfully recover the on-shell coUinear 
limit when k is coUinear with p. 

The evaluation of traces was done using the algebraic manipulation system FORM [19]. 
The expression for \M.QQ{eq* — )■ e7g)p can be easily obtained from (1) — (3) if we replace 
Pe k, Pe' Pq, k ^ Pe, Pq ^ Pe' and multiply (1) by an extra factor Cq. For the reader's 
convenience, we collect the analytic formulas for the off-shell matrix elements |A4LL(eg* 
e^q)]"^ and \M.QQ{eq* — >■ e^q)\'^ in the Appendix. 



2.2 Cross section for the prompt photon production 

According to the A;T-factorization theorem, the cross section of the process ep — )> e'yX can 
be written as a convolution of the off-shell matrix element \M.{eq* — > e^q)]"^ and unintegrated 
quark distribution fq{x, k^, //^): 

<^..,QQ(ep ^ e^X) - E / 256.3,2,^[p^^| exp(,^) l-^^^^QQ^^'i* ^ ^7g)r x 



27r 27r 27r' 



where y^'-i Vq and 0e' and 0^ are the center-of-mass rapidities and azimuthal angles of the 
outgoing electron and (anti)quark, respectively. The rapidity of the produced photon is 
given by 

^/s - me'Texp{ye') - mqTexp{yq) 



(8) 



\Pit\ 

where rrie'T and m^T are the transverse masses of the corresponding particles. We note that 
averaging the expression (7) over cj) and taking the limit k|. — )■ we recover the well-known 
LO result of coUinear parton model. We use the unintegrated quark densities in a proton 
fq{x, kl, /x^) taken in the KMR form [18]. The KMR approach is a formalism to construct the 
unintegrated parton (quark and gluon) distributions from the known conventional parton 
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distributions xa(x,/i^), where a = g or a = q. For the input, we have used the recent 
leading-order Martin-Stirhng-Thorne-Watt (MSTW) [20] parton densities. 

As it is often done for deep inelastic prompt photon production, we choose the renor- 
malization and factorization scales to be /x^ = ^Q^. In order to estimate the theoretical 
uncertainties of our calculations we vary the scale parameter ^ between 1/2 and 2 about the 
default value ^ = 1. We use the LO formula for the strong coupling constant with 
n/ = 4 active (massless) quark flavours and Aqcd = 200 MeV, such that as(M|) = 0.1232. 

The multidimensional integration in (7) has been performed by the means of Monte Carlo 
technique, using the routine VEGAS [21]. The full C++ code is available from the authors 
on requeslH. 

2.3 Fragmentation contributions 

Perturbation theory becomes nonapplicable when the wavelength of the emitted photon 
(in the emitting quark rest frame) becomes larger then the typical hadronic scale O (1 
GeV~^). Then the nonperturbative effects of hadronization or fragmentation must be taken 
into account. Accordingly, the calculated cross section can be split into two pieces 

da = (icrdircct(/W^) + (icrfragm(/i^) (9) 

with da direct ifJ''^) representing the perturbative contribution and dafj-agmifJ''^) the fragmen- 
tation contribution. In our calculations we choose the fragmentation scale fi^ to be the 
invariant mass of the quark + photon subsystem, /i^ = (p + Pi)^, and restrict (io"dircct(At^) to 
/i > M ~ 1 GeV. Under this condition, the contribution (i(Tdirect(Ai^) is free from divergences. 
We have checked that the sensitivity of our results to the choice of M is reasonably soft. As 
far as the fragmentation contribution is concerned, its size is dramatically reduced by the 
photon isolation cuts (see below). According to the estimates presented in [22], the contri- 
bution from (i(jfragm amounts to only few percent of the visible cross section. This value is 
smaller than the theoretical uncertainty in the perturbative contribution da^hect, and so, is 
neglected in our analysis. 

2.4 Photon isolation cuts 

In order to reduce huge background from the secondary photons produced by the decays 
of 7r° and 1] mesons the isolation criterion is introduced in the experimental analyses. This 
criterion is the following. A photon is isolated if the amount of hadronic transverse energy 
E^p'^, deposited inside a cone with aperture R centered around the photon direction in the 
pseudo-rapidity and azimuthal angle plane, is smaller than some value E^^^: 

ETphad ^ rpmax 
IZ/rp ^ Il/rp , 

had (10) 
(^had_^)2^(^had_^)2<^2_ 

Both the HI and ZEUS collaborations take ~ 1.0 and E^^^ ~ 1 GeV [1-3]. Isolation 
not only reduces the background from light hadron decays but also significantly reduces the 

lipatov@theory.sinp.insu.ru 
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fragmentation components connected with collinear photon radiation. It was shown that 
after applying the isolation cut (10) the contribution from the fragmentation subprocesses 
is suppressed [22]. 



3 Numerical results 

3.1 Inclusive production 

Experimental data on the deep inelastic inclusive prompt photon production come from 
both the ZEUS and HI collaborations. The differential cross sections are measured as func- 
tion of the produced photon transverse energy Ej,, pseudo-rapidity r]'^, exchanged photon 
virtuality and Bjorken variable x. The ZEUS data [1] refer to the kinematic region defined 
as 5 < < 10 GeV, -0.7 < r]^ < 0.9, > 35 GeV^, E'^ > 10 GeV, 139.8° < 6'^ < 171.8°. 
The initial electron and proton energies are Ee = 27.6 GeV and Ej, = 920 GeV. The E'^ 



and 6'^ are the energy and the polar angle of the scattered electroijj. Another, very recent 
set [2] of the ZEUS data, refers to the kinematic region 4 < < 15 GeV, -0.7 < t]^ < 0.9, 
10 < g2 < GeV^, E'^ > 10 GeV, 139.8° < 9'^ < 171.9°. The invariant mass of the 
final state hadronic system Wx > 5 GeV. The HI data [3] refer to the kinematic region 
3< E^ <10 GeV, -1.2 < t]^ < 1.8, A < < 150 GeV^, E'^ > 10 GeV, 153 < 9'^ < 177, 
y > 0.05. Additionally, the HI data on the transverse energy distributions have been divided 
into five subdivisons of i]'^, namely —1.2 < 7]"' < —0.6, —0.6 < r]'^ < 0.2, 0.2 < t]'^ < 0.9, 
0.9 < 7]"' < 1.4 and 1.4 < t]'^ < 1.8. To suppress the contribution from the elastic Compton 
scattering a cut on the mass of the final state hadronic system Wx > 50 GeV is applied. 

Our predictions on the differential cross sections as a function of transverse energy Ej^, 
pseudo-rapidity rj"', Q"^ and x variables are shown in Figs. 1 — 5 in comparison with the 
data [1-3]. Solid histograms are obtained by fixing both the factorization and normalization 
scales at the default value /i^ = Q^, whereas the upper and lower dashed histograms cor- 
respond to the scale variation as it was described above. We find that the predicted cross 
sections agree well with the HI and ZEUS data both in the normalization and shape. A slight 
underestimation of the HI data is only observed at high Q^. In general, the fc^- factorization 
predictions describe the data much better than the LO collinear calculations. 

Our calculations give for the total cross sections 47.7l};5 pb, S.ll'^oJ pb and 22.0lo:9 pb 
for the HI [3], ZEUS [1] and ZEUS [2] kinematic regions, respectively. The experimentally 
measured cross sections are 50.3 ± 1.7 (stat.) tyj (syst.) pb, 5.64 ± 0.58 (stat.) '^^'^2 (syst.) pb 
and 19.4 ± 0.7 (stat.) lJ;o (syst.) pb. The collinear LO pQCD predictions are 28.6 pb and 
5.39 pb for the HI [3] and ZEUS [1] kinematical regions. Besides that, the sole QQ contri- 
bution has been extracted in the ZEUS experiment [2] giving 12.2 ± 0.7 (stat.) t\l (syst.) pb. 
Our prediction of 14.8]'^o;5 pb is very close to this value. We find that the QQ contribution 
yelds about 65% of the total cross section. In the collinear LO calculations this mecha- 
nism gives only about 50%. The reason for the enhancement of the QQ contribution in 
the /cT-factorization approach can be seen in the fact that the unintegrated parton (quark) 
distribution effectively include a large piece of higher-order corrections. According to the 

^AU kinematic variables are given in the laboratory frame where the positive direction of the OZ axis is 
given by the proton beam. 
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observation of [3], the difference between the colhnear LO calculations and the data can 
mainly be attributed to an underestimation of the QQ contribution. 

The dependence of our results on the renormalization/factorization scale /i^ on its own 
weak, leading to about 5% uncertainty band over a wide kinematic range. We cannot 
estimate the uncertainty connected with the choice of unintegrated quark distributions (al- 
though, it may be potentially much larger) because, at the time being, they are only available 
in the KMR schem^ (see also reviews [13] for more information). 

3.2 Production in association with a jet 

To calculate the jet-associated production of prompt photons we apply the procedure used 
previously in [24]. The produced photon is accompanied by a number of partons radiated 
in the course of the parton evolution. On the average, the parton transverse momentum 
decreases from the hard interaction block towards the proton. As an approximation, we 
assume that the parton with momentum k' emitted at the last evolution step compensates 
the transverse momentum of the parton participating in the hard subprocess, i.e. k'-p ~ — k^. 
All the other emitted partons are collected together in the proton remnant, which is assumed 
to carry only a negligible transverse momentum compared to kV- This parton gives rise to 
a final hadron jet with E^r^^ = IkVI, in addition to the jet produced in the hard subprocess. 
From these hadron jets we choose the one carrying the largest transverse energy. 

Experimental data for this process have been obtained by the HI and ZEUS collabora- 
tions. The ZEUS collaboration presented the cross sections [1] measured in the kinematic 
region 5 < < 10 GeV, Ei^' > 6 GeV, -0.7 < r]^ < 0.9, -1.5 < r/j^* < 1.8, > 35 GeV^ 
E'^ > 10 GeV, 139.8° < 9'^ < 171.8° with electron energy E^ = 27.6 GeV and proton energy 
Ep = 920 GeV. The more recent HI data [3] refer to the kinematic region 3 < Ej, < 10 GeV, 
E}^' > 2.5 GeV, -1.2 < t]^ < 1.8, -1.0 < r;j^* < 2.1, A < < 150 GeV^, E'^ > 10 GeV, 
153° < 9'^ < 177°, y > 0.05 and Wx > 50 GeV with the same electron and proton energies. 

The results of our calculations are shown in Figs. 6 and 7 in comparison with the HERA 
data. One can see that the distributions measured by the HI collaboration are well re- 
produced by our calculations. However, our results overeshoot the ealier ZEUS data. A 
possible reason for that can be connected with the jet selection algorithm, being the con- 
sequence of the used approximations. Another possible reason is that our calculations are 
restricted to the parton level only, while the jet fragmentation mechanism could also be of 
importance. Concerning the colhnear NLO predictions, the shapes of the differential cross 
sections da/dEj,, da/dr]'^ and da/dQ^ are described well, but the normalization is too low, by 
about 40% [3]. The difference between the LO and NLO predictions is mostly concentrated 
at low [3]. 

The calculated total cross sections a{ep — )■ 67 + jet) = 33.2t5:o pb and 1.9io;i pb have 
to be compared with 31.6 ± 1.2 (stat.) ^'^^ (syst.) pb and 0.86 ± 0.14 (stat.) to;34 (syst.) pb, 
measured by the HI and ZEUS collaborations in the relevant kinematic regions. We find 
perfect agreement between our predictions and recent HI data. A sizeble overestimation of 
the ZEUS data is observed. It seems that the consistency of the HI and ZEUS data with each 
other can be questioned. The total cross sections calculated in the colhnear approximation 
at the LO and NLO level for HI conditions are 16.7 pb and 20.2 pb, respectively [3]. 

^Several attempts to calculate the CCFM-evolved unintegrated quark densities were made in [23]. 
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As the final point, we should mention that corrections for hadronisation and multiple 
interactions have been taken into account in the NLO analysis of the HERA data [1-3] 
performed in the framework of coUinear factorization. The correction factors are typically 0.8 
— 1.2 depending on a bin. These corrections are not taken into account in our consideration. 



4 Conclusions 

We have investigated the deep inelastic production of prompt photons at HERA in the 
/cT-factorization approach. Our study is based on the off-shell partonic amplitude eg* — ?■ e7g, 
where the photon radiation from the leptons as well as from the quarks is taken into account. 
The unintegrated quark densities in a proton are calculated using the Kimber-Martin-Ryskin 
prescription. 

We have studied both the inclusive and jet associated production of prompt photons. 
Our numerical predictions on the inclusive production cross sections arc in well agreement 
with the HI and ZEUS data. We have demonstrated that in the /cT-factorization approach 
the role of the QQ contribution is enhanced compared to the coUinear LO approximation 
of QCD. Our results for the jet associated production agree with the HI measurements but 
overshoot the ZEUS data. Probably, there exists a problem of incompatibility of these HI 
and ZEUS data with each other. 
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5 Appendix 

Here we present compact analytic expressions for the matrix elements involved in (7). 
In the massless limit, the squared off-shell matrix elements |A^LL,QQ(eg* — >■ e7g)p summed 
over the final states and averaged over the initial states are 



\M 



LL\ — 



16(A;2 - 2(pg • k)f 



7(1) 

LL 



+ 



7(2) 
LL 



+ 



2F, 



(12) 



LL 



{P^ ■ Pe'Y {P^ ■ Pef (P^ • Pe) (P^ " Pe') 



(Al) 



\Mqq\' = 



e^iPe-Pe'Y 



.(1) 



.(2) 



+ 



+ 



2F, 



(12) 



where 



= l28{p''-p,){p^-pe'){Pe-Pp) + I28{p''-Pe'){p^-Pp){p,-Pe), 



{A.2) 



{A.3) 
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F^^l = 128(p^ -pjb^ -Pejbe -Pp) + 128(p^ •Pe)(p^ -Ppjb, -PeO, (^-4) 

FlL = -64(p'^ ■ Pg)(pe • Pe'){Pe " Pp) + Q^P^ " PQ)(Pe ' Pe'){Pe " Pp)" 
64 (p'^ ■ Pe){Pq ■ pe){pc ' Pj;) - 64 (jP' • /^.^(j^q • ' Pp)- 

128(p^ • Pe)(Pg • Pe'){Pe ' Pp) + 128(p^ • Pe'){Pq ' Pe){Pe ' Pp) + 
64(p^ • Pe')(Pq ■ Pe)(Pe ' Pp) + 64(p^ • Pe'){Pq ' Pe')(Pe ' Pp)- 
64(p^ • Pp){Pq ■ Pe)(Pe ' Pe') + 64(p^ • Pp){pq ■ Pe')(Pe ' Pe')- 

128(pg ■ Pe){Pe ■ Pe'){Pe " Pp) " 128 ■ Pe')(Pe " Pe')(Pe " Pp): {^.5) 

F^'l = 128(p'^ ■Pe){Pe'-Pp) + 128(p'' -peOlPe - Pp), (^-6) 

= 128(p^ • Pe)(p^ • Pp){Pq ■ Pe') + 128(p^ • Pe')(P^ " Pp) (P<? " Pe) + 
128(p'^ • A;)(pg ■Pe)(Pe' -Pp) + 128(p'^ • k){pq ■ Pe'){Pe ■ Pp)- 
l2S{p^ ■ Pp){Pq ■ Pe){Pe' ■ k) - 128 (p"^ ■ Pp){Pq ■ Pe')(Pe " k)- 

6A{pq ■ Pe){Pe' ' Pp)k'^ " 64(^5 • Pe')(Pe " Pp)^^ (^-7) 

FqQ = 64(p^ • Pq){pq ■ Pe){Pe' ' Pp) + 64(p^ • Pq){pq ■ Pe'){Pe ' Pp) + 
64(p^ • Pq){Pe ■ k){pe> ■ Pp) + 64(p^ • Pq){pe ■ Pp){Pe' ' k) + 
64(p^ • Pe)(Pq ■ Pe'){Pq ' Pp) " 64(p^ • Pe){Pq ' k)(pe' ■ Pp) + 
64(p^ • Pe')(P<Z • Pe){Pq ' Pp) " 64(p^ • Pe'){Pq ' A;)(pe ' Pp)" 
QA{p"' ■ k){pg ■ Pp){pe ■ Pe') - 128(p^' ■ Pp){pq ■ Pe){Pq " Pe')" 
64(p'^ ■ Pp)(Pg ■ Pe)(Pe' ■ k) - 64(p'^ ■ Pp){Pq ■ Pe'){Pe " ^) + 
U{p^ ■ Pj,){p„ ■ k){pe ■ Pe') - 64(pg ■ Pe){Pq ■ k){pe' ■ Pp)- 
64(Pg • Pe){Pq ■ Pp){Pe' ' k) - 64(pg • Pe')(Pg " k){pe ■ Pp)- 

64(Pg • Pe'){Pq ■ Pp)(Pe ' k). (A.8) 
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Figure 1: Differential cross sections of the inclusive deep inelastic prompt photon production 
as a function of E^, 7]^ and Q'^ calculated at 3 < £^5^ < 10 GeV, -1.2 < r^^ < 1.8, A < < 
150 GeV^, E'^ > 10 GeV, 153° < 9'^ < 177°, y > 0.05 and Wx > 50 GeV. The sohd histogram 
corresponds to the default scale ji'^ — Q'^, whereas the upper and lower dashed histograms 
correspond to scale variations described in the text. The experimental data are from HI [3]. 
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Figure 2: Differential cross sections of tlie deep inelastic inclusive prompt photon production 
as a function of e:^ and 77^ calculated at 3 < < 10 GeV, -1.2 < 77T < 1.8, 40 < < 
150 GeV^, E'^ > 10 GeV, 153° < 9'^ < 177°, y > 0.05 and Wx > 50 GeV. Notation of the 
histograms is as in Fig. 1. The experimental data are from HI [3]. 




Figure 3: Differential cross sections of the deep inelastic inclusive prompt photon production 
as a function of E]. and r]^ calculated at 5 < < 10 GeV, -0.7 < r]^ < 0.9, > 35 GeV^, 
E'^ > 10 GeV and 139.8° < 6'^ < 171.8°. Notation of the histograms is as in Fig. 1. The 
experimental data are from ZEUS [1] . 
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Figure 4: Differential cross sections of the inclusive deep inelastic prompt photon production 
at HERA as a function of E^, rp , and x calculated at 4 < i?^ < 15 GeV, -0.7 <rp < 0.9, 
10 < g2 < 350 GeV^, E'^ > 10 GcV, 139.8° < d'^ < 171.9° and Wx > 5 GeV. Notation of 
the histograms is as in Fig. 1. The experimental data are from ZEUS [2]. 
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Figure 5: Differential cross sections of the deep inelastic inclusive prompt photon production 

as a function of calculated for different subdivisons of the r)'^ range. Notation of the 
histograms is as in Fig. 1. The experimental data are from HI [3]. 
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Figure 6: Differential cross sections of tlie deep inelastic prompt photon and jet associated 
production as a function of E^, rj^ and calculated at 3 < < 10 GeV, E^r^' > 2.5 GeV, 
-1.2 < < 1.8, -1.0 < Tyj^t < 2.1, 4 < < 150 GeV^, E'^ > 10 GeV, 153° < 9'^ < 177°, 
y > 0.05 and Wx > 50 GeV. Notation of the histograms is as in Fig. 1. The experimental 
data are from HI [3]. 
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Figure 7: Differential cross sections of the deep inelastic prompt photon and jet associated 

production as a function of Ej,, rj'^, E^r^* and r/j*^* calculated at 5 < E^ < 10 GcV, i?;^* > 
6 GeV, -0.7 <r]^ < 0.9, -1.5 < < 1.8, > 35 GcV^, E'^ > 10 GeV and 139.8° < 9'^ < 
171.8°. Notation of the histograms is as in Fig. 1. The experimental data are from ZEUS [1]. 
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